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Gamma-ray emission from novae related to positron annihilation: constraints 
on its observability posed by new experimental nuclear data 
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ABSTRACT 

Classical novae emit 7-ray radiation at 511 keV and below, with a cut-off at around 
(20-30) keV, related to positron annihilation and its Comptonization in the expanding 
envelope. This emission has been elusive up to now, because it occurs at epochs well 
before the maximum in optical luminosity, but it could be detected by some sensitive 
intrument on board a satellite, provided that the nova is close enough and that it is 
observed at the right moment. The detection of this emission, which is a challenge for 



the now available and for the future 7-ray instruments, would shed light into the physi- 
cal processes occurring in the early phases of the explosion, which are invisible in other 
lower energy ranges. A good prediction of the emitted fluxes and of the corresponding 
detectability distances with different instruments relies critically on a good knowledge 
of reaction rates relevant to ^^F destruction, which have been subject to a strong revi- 
sion after recent nuclear spectroscopy measurements. With respect to previous results, 
smaller ejected masses of ^^F are predicted, leading to smaller emitted fluxes in the 
(20-511) keV range and shorter detectability distances. 

Subject headings: gamma-rays: observations — novae, cataclysmic variables — nuclear 
reactions, nucleosynthesis, abundances 



1. Introduction 

The role of novae as potential gamma-ray emitters has been mentioned long ago (Clayton 
& Hoyle 1974, Clayton 1981, Leising & Clayton 1987), and analyzed in detail on the basis of 
theoretical models just recently (Gomez-Gomar et al. 1998). In many of the historical studies, 
more attention has been paid to the long lasting emission, related to the decay of some medium- 
lived isotopes, such as '''Be (r=77days), ^^Na (r=3.75yr), and to the possible contribution of novae 
to the Galactic content of ^^Al (r = 1.04 x 10^ yr). Up to now, no positive detection neither of the 
478 keV emission (Harris et al. 1991), related to ^Be decay, nor of the 1275 keV one (lyudin et al. 
1995), related to ^^Na decay, has been obtained. The 1809 keV line associated with the galactic ^^Al 
was detected some years ago by the HEA03, High Energy Astrophysics Observatory (Mahoney et 
al. 1982), but observations made with the COMPTEL (Compton Telescope) instrument on board 
the CGRO satellite (Compton Gamma-Ray Observatory) seem to indicate that the ^^Al emission is 
more related to a young population of massive stars, although a nova contribution cannot be ruled 
out yet (Diehl et al. 1995, Prantzos & Diehl, 1996). The most powerful emission in gamma-rays 
from classical novae is not the one related to the 478, 1275 or 1809 keV lines, but that at 511 keV 
and below (down to ~ (20-30) keV), which originates from electron-positron annihilation (Hernanz 
et al. 1997a, b, Gomez-Gomar et al. 1998), and consists of a line at 511 keV and a continuum 
below it. The 511 keV line comes from the direct annihilation of positrons and from the positronium 
(in singlet state) emission, whereas the continuum comes from both the positronium continuum 
(triplet state positronium) and the Comptonization of photons emitted in the 511 keV line. The 
main contributors to positrons in nova envelopes are the short-lived ^^N (r=862s) and ^^F (r=158 
min), besides the long-lived ^^Na (this one only in oxygen- neon novae). The problem with this 
emission is its very short duration, making its detection very difficult. Furthermore, because of the 
short lifetimes of ^^N and ^^F , its maximum happens well before the maximum in visual light, 
which still enhances the difficulty of detection. When ^^N decays, the nova envelope is still too 
opaque; therefore, ^^F plays the most important role, since it decays appreciably when the envelope 
begins to be transparent enough to allow for the gamma-ray photons to be transported through 



the whole envelope, and then be emitted into space. 

Several attempts are being carried out to detect the annihilation gamma-ray emission from no- 
vae, with instruments like BATSE (Burst and Transient Source Experiment) on board the CGRO, 
or TGRS (Transient Gamma Ray Spectrometer), on board the WIND satellite. SPI, the future 
spectrometrer for INTEGRAL (International Gamma- Ray Laboratory), could in favorable condi- 
tions detect that emission. The detectability distances for these instruments, or for any future ones, 
rely on the detailed properties of novae expanding envelopes and, in particular, on their ^^F con- 
tent, which strongly depends on the rates of the nuclear reactions involved in its synthesis. A good 
knowledge of the gamma-ray output of novae is needed to elucidate the possibilities of detection 
and to put constraints on the theoretical models from the observations. 



2. Gamma-ray emission of novae and ^^F synthesis 

The flux emitted by a classical nova in the gamma-ray range below 511 keV depends on the 
amount of ^^F at the epoch where the envelope begins to be transparent to gamma-rays. Therefore, 
profiles of ^^F along the expanding envelope at different times are needed to compute the spectrum 
at different epochs after the outburst (or in an equivalent way, the light curves for different gamma- 
ray bands at and below 511 keV). Of course, the whole spectrum depends also on the abundances of 
the other isotopes and on the density and velocity profiles; these quantities determine the number of 
interactions of photons with the ejecta through Compton scattering, photoelectric absorption and 
pair production. Theoretical models of the gamma-ray emission from individual classical novae, 
based on the nucleosynthetic yields obtained from a hydrodynamical code that follows the explosive 
phase and the related nucleosynthesis (Jose & Hernanz 1998), predict an intense line emission at 
511 keV plus a continuum below (Gomez-Gomar et al. 1998). The emission appears between 5 and 
6 hours after maximum temperature, i.e., some days before the maximum in visual luminosity. The 
continuum displays always a cut-off at low energies, which is related to photoelectric absorption, 
which has larger cross sections than Compton scattering at low energy, acting as a sink of the 
Comptonized photons. The cut-off is located at ~20 keV, for the carbon-oxygen (CO) and ~30 
keV for the oxygen-neon (ONe) rich envelopes typical of nova models. 

The main nuclear path leading to ^^F synthesis in novae belongs to the hot CNO (carbon- 
nitrogen-oxygen) cycles. Since in both CO and ONe novae the initial abundance of ^^O is large, 
^^O is the main source for ^^F formation, through two possible chains of reactions: ^^O (p,7)^^F 
(p,7)^^Ne(/3+)^^F and ^"^O {p,^y'^F (/3+)^'^0(p,7)^^F . However, ^^F yields are severely constrained 
by its destruction mode, whatever the production channel is. During the thermonuclear runaway, 
^^F destruction by beta decays can be neglected when compared to its destruction by proton 
captures. Because of the low alpha emission threshold, ^^F (p,a)^^0 is faster than ^^F (p,7)^^Ne, 
and hence it is the main destruction channel of ^^F (see Figure 1). The rates for ^^F (p,7) and 
""^^F (p,a) generally adopted in nova models are those from Wiescher & Kettner (1982), which were 
estimated at a time where the relevant experimental data were scarce (i.e., very limited spectroscopic 



data was available for -^^Ne). Thanks to recent measurements, in particular to those of Graulich 
et al. (1997), Utku et al. (1998), and references therein, it is possible now to use updated rates. 
Two resonance strengths have been measured directly (Coszach et al. 1995, Graulich et al. 1997) 
in ^^F (p,a) using a radioactive ^^F beam. Of these two, the one lying at higher energy is so 
broad that its tail alone gives a contribution larger than the Wiescher &: Kettner (1982) rate at 
nova temperatures. Nuclear spectroscopic studies performed by Utku et al. (1998) have led to 
the localisation of several new ^^Ne levels to which should correspond new resonances in ^^F +p. 
They have calculated updated ^^F (p,7) and ^^F (p,a) rates including the two directly measured 
resonances plus the contribution of the expected ones. To calculate this last contribution, they used 
information about ^^F analog levels supplemented with estimated quantities. A comparison of the 
new rates with those of Wiescher and Kettner is shown in Figure 1; in the range of temperatures 
relevant to nova explosions, the new rates are quite larger than the older ones. Although the 
uncertainties associated with the rates of ^^F destruction, ^^F (p,7) and ^^F (p,a), are far from 
being settled in the nova temperature range, it seems clear that the rates that had been used up 
to now (Wiescher & Kettner, 1982) were underestimated. 

In Table 1 we show the most relevant properties of some nova models computed by means 
of the hydrodynamic code described in Jose & Hernanz (1998), with the nuclear reaction network 
described in Jose et al. (1999) (recommended rates) and the new Utku's rates for ^^F (p,7) and ^^F 
(p,a) reactions. We have chosen models representative of the two nova types, CO and ONe, and 
with masses 1.15 Mq in the CO case and 1.15 and 1.25 Mq in the ONe case, in order to analyze 
the influence of the composition and mass of the underlying white dwarf on ^^F synthesis. The 
accretion rate in all the cases is 2 x 10^^*^ MQ.yr"^ and the degree of mixing with the core is 50% 
(see Jose &: Hernanz 1998 for a discussion of the initial mixing effect). Similar ejected masses of ^^N 
and ^^F are obtained in all the models, which is not the case for other radioactive nuclei relevant to 
7-ray emission with longer timescales, such as ^Be (mainly produced in CO novae) and ^^Na and 
^^Al (mainly synthesized in ONe novae). When compared with results with the old rates for ^^F 
+p reactions, a general reduction, by a factor between 10 and 20, of ^^F yields is present, whereas 
the maximum temperature and the mean kinetic energy are almost unchanged, since the ^^F +p 
reactions do not affect the energetics of the nova explosion. Concerning the general nucleosynthesis, 
the changes are not very relevant, except for ^^F, which is reduced typically by a factor of 10. 

The gamma-ray spectrum 12 hours after peak temperature is shown in Figure 2, for a CO 
nova of 1.15 Mq and an ONe nova of 1.25 Mq at a distance of 1 kpc. The 511 keV line from our 
models is slightly blueshifted (with a maximum energy of 517 keV) during the first phases, and has 
a moderate width (from ^^5 to ~8 keV, FWHM, depending on the model and on the epoch). Light 
curves in the 20-511 keV band and in the 511 keV line alone have been computed. A summary of 
their fluxes at different epochs after peak temperature, for a distance of 1 kpc, is shown in Table 
2. There is an early maximum at around 1 hour, which is related to ^^N decay. This maximum 
has very short duration and is very dependent on the distribution of ^^N in the outer layers of the 
envelope (which are the only ones seen at these early epochs). A larger abundance of ^^N in the 



outer layers of the CO nova with respect to the ONe one leads to the larger flux in the CO case 
(see Table 2). The flux reaches its second and broader maximum, related to ^^F decay, at around 
6 hours after peak temperature. It is important to notice that the flux emitted in the continuum, 
with a maximum around 50 keV, is larger than that in the line at all times. Also, the clear cut- 
off at ~20 keV, together with the short duration of the whole emission and its early appearence 
before maximum visual luminosity, rules out the interpretation of the early hard X-ray emission 
observations of Nova Herculis 1991 (Lloyd et al. 1992), and of the recent Nova Velorum 1999 (Orio 
et al. 1999, Mukai &: Ishida 1999), as a Comptonization of the gamma-ray emission (Livio et al. 
1992, Starrfield et al. 1992). For the purpose of comparison, the gamma-ray spectrum of the CO 
nova model, but computed with the old ^^F -|-p rates, is also shown in Figure 2. A reduction by a 
factor of ~10 at all the energies in the relevant range is obtained in the new models, since the fluxes 
are smaller than the previous ones by the same factors than ^^F ejected masses (between 10 and 
20). This difference appears during the first 24 hours, when ^^F is the main source of positrons; 
but later on, when only ^^Na remains as a positron source, the new and the old fluxes are quite 
similar, being much smaller in CO novae because they are almost devoid of ^^Na . 



3. Discussion and conclusions 

The potential detection of the electron-positron annihilation emission from classical novae is 
faced with various challenging problems: besides of the relatively small fluxes, the short duration 
and the early appearence of the emission make its detection quite hard. Large field of view in- 
truments, like BATSE, are the more appropriate to detect this type of emission. Its capability to 
observe all the sky, together with its high sensitivity in the low-energy range make BATSE an ideal 
instrument to detect this emission, although the new fluxes presented here restrict the detectability 
distances to less than 3 kpc. There is an ongoing project of analysis of BATSE background data 
for nearby novae that have exploded since CGRO launch. 

Another instrument which could be able to detect the annihilation gamma-ray emission from 
novae, because of both its wide field of view and its high spectral resolution, is the Transient 
Gamma Ray Spectrometer (TGRS) on board the WIND satellite (Teegarden et al. 1996). Very 
recently, Harris et al. (1999) have tried to detect the annihilation line emission at 511 keV from 
the five known Galactic novae in the instrument field of view, during the time interval 1995-1997. 
No definite detections were made, but the method applied was shown to be sensitive enough to 
detect novae out to about 2.5 kpc. Since the conclusions of that paper were based on models in 
Gomez-Gomar et al. (1998), they should be revised and the corresponding distances should be 
reduced by a factor ~ -v/IO- 

Finally, we should mention the future spectrometer for INTEGRAL, SPI, which could be able 
to detect the 511 keV line up to 3 kpc with the fluxes we obtain now. A similar detectability 
distance is predicted for the continuum. It is worth noticing that both TGRS and SPI have very 
good spectral resolution, because of their germanium detectors; this of course is favorable for the 



detection of the 511 keV line (although its non negligible width worsens the sensitivity appreciably) 
but not so interesting for the continuum, where novae emit a larger flux. 

The detection by some present or future instrument of the positron annihilation emission from 
novae would help to understand the mechanism of the explosion and the dynamics of the early 
expansion stage. 

This research has been partially supported by the CICYT-P.N.I.E. (ESP98-1348), by the DG- 
ICYT (PB97-0983-C03-02; PB97-0983-C03-03), by the CIRIT (GRQ94-8001) and by the PICS 
319. 
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Fig. 1. — Comparison between the Utku et al. (1998) rates and those from Wiescher & Kettner 
(1982), as well as comparison between -"^^F (p,a) and ^^F (p,7) rates, for a wide temperature range. 

Fig. 2. — Gamma-ray spectrum at t=12 h after peak temperature, for a CO nova of 1.15 Mq (solid 
line) and an ONe nova (dashed line). The old spectrum from Gomez-Gomar et al. (1998) for the 
CO nova model is also shown for comparison (dotted line). 



Table 1: Main properties of the ejecta 1 hour after peak temperature. Tmax is the peak temperature, 
< -E'kin > the mean kinetic energy, in units of 10^^ ergs. The ejected masses of ^^N and ^^F are in 







Nova type M„d (Mq ) T^^^jK) < ^kin > Mej(^^N ) Mej(^^F ) 

CO 1.15 2.05 X 10^ 1.1 2.3 x 10"^ 2.6 x lO^^ 

ONe 1.15 2.31 x 10^ 1.5 2.9 x 10"^ 5.9 x 10"^ 

ONe 1.25 2.51 x 10^ 1.5 3.8 x 10"® 4.5 x 10"^ 



Table 2: Fluxes, in photons. s~^.cm~^, at different times after peak temperature, in two energy 
bands in the continuum and in the 511 keV line (FWHM=8 keV), for a CO nova of 1.15 Mq and 
an ONe nova of 1.25 M.q , at a distance of 1 kpc. The fraction of escaping energy in gamma-ray 
photons is displayed in the last column. 

t (hours) Model F(2o_250)keV F(25o-5ii)kcV Fsiikcv line fraction escaping (in %) 

1 CO 2.5 X 10"^ 9.6 X 10^2 2.2 x 10"^ 0.4 

ONe 2.2 X 10^2 93 ^ 10"3 ^_g ^ ^^O^^ 0.02 

3 CO 5.3 X 10^3 1.8 X 10^3 3.6 x 10"^ L9 

ONe 1.6 X 10^2 7 2 x 10"^ 1.6 x lO'^ 4.2 

6 CO 8.1 X 10-3 2.9 X 10-3 6.0 x 10"^ 10 

ONe 1.7 X 10-2 8.8 x 10-^ 1.9 x 10-^ 17 

12 CO 2.5 X 10-3 1.2 X 10-3 3.0 x 10-"^ 39 

ONe 3.8 X 10-3 2.6 x 10-3 6.7 x 10-^ 47 

24 CO 2.5 X 10-^ 3.0 x 10-^ 8.4 x 10-^ 77 

ONe 5.2 X 10-^ 7.4 x 10-^ 2.0 x 10-^ 80 

48 CO 1.5 X 10-"^ 2.3 x 10-^ 2.0 x 10-^ 95 

ONe 1.1 X 10-^ 2.7 x 10-^ 8.2 x 10-^ 94 
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